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Abstract. In this second paper, using N = 3 polarized electrons (spinless fermions) interacting via a U/r 
Coulomb repulsion on a two dimensional L x L square lattice with periodic boundary conditions and nearest 
neighbor hopping t, we show that a single unpaired fermion can co-exist with a correlated two particle 
Wigner molecule for intermediate values of the Coulomb energy to kinetic energy ratio r a = UL/(2t\/ tvN). 
This supports in an ultimate mesoscopic limit a possibility proposed by Andreev and Lifshitz for the 
thermodynamic limit: a quantum crystal may have delocalized defects without melting, the number of sites 
of the crystalline array being smaller than the total number of particles. When L = 6, the ground state 
exhibits four regimes as r s increases: a Hartree-Fock regime, a first supersolid regime where a correlated 
pair co-exists with a third fully delocalized particle, a second supersolid regime where the third particle is 
partly delocalized, and eventually a correlated lattice regime. 



PACS. 71.10.-w 
Electron solids 



Theories and models of many-electron systems - 73. 21. La Quantum dots - 73.20.Qt 



1 Introduction 

In 1969, it was conjectured by Andreev and Lifshitz 
that at zero temperature, delocalized defects may exist 
in a quantum solid, as a result of which the number of 
sites of an ideal crystal lattice may not coincide with the 
total number of particles. Originally, this conjecture was 
proposed for three dimensional quantum solids made of 
atoms (He 3 , He 4 , . . .) which do not interact via Coulomb 
repulsion. We re-visit such a possibility for electron solids 
with long range Coulomb repulsion in two dimensions. 
The motivation to re- visit nowadays this issue comes from 
questions raised by the physics of electrons in Si MOS- 
FETs and similar 2c? field effect devices. An unexpected 
low temperature metallic behavior [5] has been observed at 
intermediate values of the Coulomb energy to kinetic en- 
ergy ratio r s , which remains unexplained. Another actual 
motivation is given by the promising perspectives opened 
by trapped cold ion systems, where one can study how a 
Wigner molecule becomes a quantum fluid when the ions 
are squeezed 0- 

In a first paper 0|, the supersolid phase conjectured 
by Andreev and Lifshitz was introduced, together with 
a related variational approach using a fixed number of 
fermions BCS wave function j^j of Bouchaud et al. The 
question is to know if a system of unpaired electrons with 
a reduced Fermi energy can co-exist with an ordered ar- 
ray of charges, the number of sites of the crystalline array 
being smaller than the total number of electrons. In Rcf. 



0], this question was investigated using N = 4 spinless 
fermions interacting via a U/r Coulomb repulsion in a two 
dimensional 6x6 square lattice with periodic boundary 
conditions (BCs) and nearest neighbor hopping t. It was 
observed that for intermediate ratios r s — U L / (2t^J it N) 
(typically for 10 < r s < 28), the ground state (GS) is 
in a mixed state, where unpaired delocalized fermions co- 
exist with a strongly paired, nearly solid assembly. From 
the study of the different inter-particle spacings as r s in- 
creases, it was concluded that before having full Wigner 
crystallization, a floppy three particle Wigner molecule is 
formed, while the fourth particle remains delocalized. We 
consider in this second paper the ultimate limit N = 3, 
where it is still possible to exactly study if a correlated 
two particle molecule can co-exist with a third delocalized 
particle. As in the first paper, the study is restricted to 
fully polarized electrons (spinless fermions) having anti- 
symmetric orbital wave functions. A study involving the 
spin degrees of freedom can be found in Ref. [B]. 

The paper is organized as follows. Once the lattice 
model is defined in Sec. [21 the three regimes characterizing 
the formation of a two particle Wigner molecule (2PWM) 
on an empty periodic lattice are summarized in Sec. 
and in App. fA] a weak coupling Fermi regime, a corre- 
lated Wigner regime with harmonic oscillatory motions of 
the particles around equilibrium, and a correlated lattice 
regime where these oscillatory motions become damped by 
the lattice. The weak coupling Fermi limit and the strong 
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coupling correlated lattice limit of the three particle sys- 
tem are described in Sec. 01 and Sec.03 An additional dis- 
cussion of the correlated lattice limit when L — ► oo is given 
in App.[5]both for the zero density limit (keeping N = 3) 
and for the constant density limit (taking N = L 2 /9 — 1). 
The threshold rf above which a weak coupling expan- 
sion in powers of U/t and the threshold r* under which a 
strong coupling lattice expansion in powers of t/U cease 
to be valid are defined in Sec. H3 Taking L = 6, one gets 
the range of intermediate values rf f=a 6 < r s < r* « 180 
where the GS structure is non-trivial. A simple ansatz 
for the GS wave function, first introduced in Ref. 0, is 
studied in Sec. corresponding to a 2PWM co-existing 
with a third particle which remains partly delocalized in 
the direction parallel to the 2PWM. The ansatz combines 
two possible directions for the 2PWM and a delocalized 
center of mass. This defines the concept of a partially 
melted Wigner molecule (PMWM) near the lattice limit, 
and describes the three particle GS at r s s» 40, when the 
2PWM oscillations are taken into account using a lattice 
t/U expansion. In Sec. 0we show that when r s is fur- 
ther decreased, the GS is made of a floppy 2PWM with 
large oscillatory motions which are not damped by the lat- 
tice, co-existing with a third fully delocalized particle. The 
unpaired particle simply provides a uniform background 
density for the 2PWM. The total momentum K2 of the 
2PWM and the momentum k3 of the third particle satisfy 
the conservation of the total momentum K = K2 + k3. 
This yields different possible combinations which contain 
more than more than 90 % of the exact GS when r s ss 10. 
In Sec. we summarize the four regimes found for the 
three particle GS when L = 6, pointing out the role of 
the lattice and raising the main question: Are those ob- 
served supersolid GSs the ultimate mesoscopic trace of 
a thermodynamic supersolid phase proposed by Andreev 
and Lifshitz, consisting of a N' electron solid co-existing 
with a N — N' electron fluid, out of a total number N 
of electrons? The GS nodal structure and the occupation 
numbers Pk on the reciprocal lattice are studied in App. 
Oand in App.EI 



2 Lattice model 

The Hamiltonian of the L x L square lattice model with 
periodic BCs is the same as in Ref. 0]. Denoting cj, c- } 
the creation, annihilation operators of a spinless fermion 
at the site j, it reads: 



H 



(j.j'> 



u 



E 



(1) 



The effective mass being to*, t — h 2 /(2m*a 2 ) is the hop- 
ping term between nearest neighbors, and U — e 2 /(ea) is 
the Coulomb interaction between two fermions separated 
by a lattice spacing a in a medium of dielectric constant 
e. The Coulomb energy to Fermi energy ratio r s becomes 
in this lattice model 

(2) 



2U 



where v = N/L 2 . For L = 6, this gives r s — 0.98U/t if 
N = 3 and r s = 1.20U/t if N = 2. 

Without disorder, the Hamiltonian is more conve- 
niently written using the operators (du) creating (anni- 
hilating) a spinless fermion in a single particle plane wave 
state of momentum k. The Hamiltonian becomes: 

H = —2t (cos k x + cos k y ) d^dk 

k 

+U J2 nq)4 +q 4'-<A'dk (3) 

k,k',q 



where 



T w v 1 cosc 3j 



'jo 



(4) 



The distance djo is defined as the shortest distance be- 
tween the sites j and of the square lattice with periodic 
BCs: 

d i0 = ^1 min(j x ,L - j x ) 2 + mm( j y ,L - j y ) 2 . (5) 

The states of different total momenta K are decoupled. 
Moreover, since the Coulomb repulsion is a two-body in- 
teraction, only states of same K having N—2 k in common 
out of N are directly coupled. When N > 3, this means 
that the Hamiltonian matrix of a subspace of given K is 
sparse. 



3 Formation of a two-particle Wigner 
molecule on an empty square lattice 

Before studying the three particle problem, we summarize 
the three regimes characterizing the two particle problem 
on an empty lattice. Firstly, it allows us to introduce the 
value r* above which the lattice effects become important. 
Secondly, it will be useful for analyzing in Sec.0the three 
particle GS in terms of a two particle Wigner molecule 
created on the uniform background provided by a third 
delocalized particle. 

Following Ref. |S] , we consider the relative fluctuations 



V(r 2 ) ~ W 
(r) 



(6) 



of the distance r between the two particles. This gives 
three regimes: 

— For r s < rf = 7r 3 / 2 , the fluctuation u r keeps essen- 
tially its non-interacting value, up to some negligible 
perturbative corrections. 

— For rf < r s < r* cx L 3 , u r decays as rj a , the two par- 
ticles beginning to form a correlated Wigner molecule, 
with an oscillatory motion of the particles around the 
equilibrium position of the molecule. The interaction 
Urijny /2dyy having a cusp at the equilibrium position 
if one defines dyy as previously, the oscillations are not 
harmonic and one gets a ~ 0.31 8 . Smearing this 
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Fig. 1. Relative fluctuation u r of the inter-particle distance r 
for two particles on an empty L x L square lattice (minus a 
finite size correction 0.9/(L+ l) 2 ) as a function of r s for L — 6 
(dashed line) and L = 12 (continuous line). The smeared inter- 
action (Ea. l49H has been taken. The dotted line corresponds to 
the u r oc r^ 1 ^ 4 behavior derived in App.fA] Notice the break- 
down of the universal scaling when r s > r*(L). 



Fig. 2. The GS energies E for K = (0,0) (thick line) and 
K = 2tt/6(1, 1) (thick dashed line) as a function of r s . The thin 



dotted lines give the behaviors E. 



K=(0,0) 



-10*+ 1.2933(7 and 



Sk1 2 W6(i i) = — + 1.249717 obtained at the leading order 
of a (7/t expansions (Eas. ITTl and ll2|l . 



4 The Fermi limit 



cusp, one recovers harmonic oscillations, and a = 1/4, 
as shown in App. [3 This is the harmonic regime first 
considered by Wigner 0| (see also Ref. ^DJ) I0r the 
continuum electron gas. 

If one subtracts a finite size correction 0.9/(1/ + l) 2 
(L even) from u r , one obtains a universal scaling law 
which depends only on r s w 0.2UL/t (N — 2) with 
a crossover from independent particle motion towards 
correlated motion at rf . 

When r s reaches a higher threshold r* , the oscillations 
of the inter-particle spacing become of the order of 
the lattice spacing a and a lattice expansion in powers 
oit/U becomes valid. The continuum-lattice crossover 
occurs when Ar « 2/3 (in units of a). At/U expansion 



giving: 



Ar = -L, 



(L-l) 



2tt 

T 



(7) 



one gets a lattice threshold r* oc L 3 . For N = 2 and 
L = 6, r* w 100. Above r*, u r oc ti/C7 is no longer a 
universal function of the ratio r s oc UL/t. 



The three behaviors of u r are given in Fig. 4 of Ref. jS] 
for various values of L if one takes the Coulomb interaction 
which we assume in this work. The relative fluctuations u r 
yielded by a smeared Coulomb interaction (Eg .1151 of App. 
lAl are shown in Fig. H] for L = 6 and 12. One can see that 
when r s < r*, u r — 0.9/(L + l) 2 is a universal function 



of r s with a Fermi- Wigner crossover at r s 



r 3 / 2 , and 



that this universal regime ceases to be valid due to lattice 
effects when r s exceeds r* 



We begin to study three particles on 6 x 6 periodic lattice 
when U — > 0. In this limit, the eigenstates \\P(r s = 0)) are 
plane-wave states: 



l*(r.=0)) =ct iC l 2 cl 3 \0) 



(8) 



|0) being the vacuum state. The GS energy —lOt has a 
sixfold degeneracy. A basis of this degenerate eigenspace 
can be built using two states of total momentum K = 
(0,0), given by 



r t t 



(i,o) c V(-i,o) |0> 



(9) 



and its x «-> y-symmetric counterpart, and four states 



c V(±i.o) c ap(o,±i) 1°) 



(10) 



of total momenta K = 27r/6(±l, ±1). 

When U/t is small, one can use perturbation theory to 
determine which of those six states have the lowest energy 
when one switches on U. At first order, the corrections 
AE^ to the GS energy — lOt are given by the diagonal 
elements of the interaction matrix (the two K = (0, 0) 
states being decoupled due to the additional x ^> y sym- 
metry). One gets for K = (0, 0) 



2tt 



AE (I) 

_^=6F(0,0)-4F^,0 
and for K = (2tt/6, 2tt/6): 



(11) 



AE- 



K 



U 



(,170.0)- XV ( ^,0 
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where the V(q x , q y ) are given by Eq. 01 

The sixfold degeneracy is partly removed by U, the 
four states of total momenta K = 2tt/6(±1,±1) having 
a smaller energy than the two states of K = (0,0). One 
can compare in Fig. the exact behaviors and the first 
order perturbative expansions fEas. ITTI and I12|l for K = 
2tt/6(±1, ±1) and K = (0, 0) total momenta. 



5 The correlated lattice limit 

We now consider the limit t — > 0. Without hopping t, the 
three particles stay localized on three different lattice sites, 
forming configurations which can be ordered by increas- 
ing Coulomb energy. For the low energy configurations, 
the inter-particle spacings are as large as it is possible on 
a periodic square lattice. The first configurations of min- 
imum Coulomb energy are given in Fig. [3] Without dis- 
order, the L 2 sites are equivalent, and L 2 = 36 identical 
configurations can be put on the L x L lattice unless an 
extra symmetry or the periodic BCs reduces this number. 
This yields large degeneracies when t = 0. For instance the 
GS degeneracy is equal to 36, the states being triangles 
of Coulomb energy 0.9023C/, having different locations or 
orientations (see Fig. EJ. 

This large degeneracy can be partly broken by a hop- 
ping term t ^ 0. This can be studied using a perturbation 
theory starting from the t = triangles and taking t as 
a perturbation. The first correction to Coulomb energy of 
the L 2 triangles is given at the second order. One gets 
a uniform shift cx t 2 /U which does not remove the L 2 
degeneracy. For the 36 triangles (f3 = 1 . . . 36) of energy 

£^ 0) = 0.9023C/, one gets the same second order correc- 
tion: 



E 



(<%iffii<^><<^iffii^ 



0r</ 



(4 0) 



E a ) 



(13) 



At the third order, two processes become possible for 
N = 3 and L — 6: either the N — 3 particles hop to- 
gether by one lattice spacing in the same direction, such 
that the center of mass of the corresponding triangle is 

0.9023 U 0.9269 U 0.9860 U 1.0407 U 




1.0519 U 1.0547 U 1.1106 U 



1.1324 U 



Fig. 3. Some of the low energy configurations with their elec- 
trostatic energy on L = 6 lattice. 
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Fig. 4. The L 2 — 36 triangles of Fig.l^can be put on an effec- 
tive L x L periodic lattice. The four sectors correspond to the 
four possible orientations of the triangles (shown in the corner), 
each oriented triangle having L/2 x L/2 possible locations on 
the original L x L periodic lattice. The first neighbor hops h 
and third neighbor hops r correspond to a t/U expansion at 
the order N = L/2 = 3. 

translated by the same hop (hopping term h cx t N /U 1 ^^ 1 ) 
or one particle hops over a scale L/2 = 3 (hopping term 
r tx t L l 2 /U L ^ 2 ^ 1 . This L/2 = 3 hop couples two triangles 
having in common two sites, i.e changes the orientation 
of the triangle. Those two processes are possible at the 
same order when A" = L/2, which is our case. Inside the 
subspace spanned by the 36 triangles, the translational 
invariance is recovered at the order N = L/2 = 3. Each 
eigenstate can be now labeled by its quantized total mo- 
mentum K. This partly removes the 36 degeneracy of the 
triangles. The matrix elements of the 36 x 36 secular ma- 
trix are given at the order A^ = L/2 = 3 by: 



M 



(in) 



0,0' 



E 



(% \H X | V a ) (V a | V a .) (V a |ff! | 



(E®-E a )W-E, 



(0) 



(14) 

where (3, (3' labels two different triangles. The diagonal- 
ization of this matrix is easy, since we can order the 36 
triangles of a 6 x 6 lattice as indicated in Fig. 0] A site j of 
this effective lattice corresponds to a triangle, and Dj (Dj) 
are the corresponding creation (annihilation) operators 



c\c). 



Jy+3) (jx+3,jy)' 



(15) 



in) 



One gets the effective Hamiltonian H e ^ 



H 



eff 



(JJ'> 



<jj'>; 



(16) 

which is identical to a one particle Hamiltonian describing 
the motion of a single particle on a 6 x 6 square lattice with 
periodic BCs, with first neighbor hopping matrix element 
h and third neighbor hopping matrix element r. h and r 
are given by the corresponding matrix elements of . 

For A^ = 3 and L — 6, the eigenenergies of are: 



AE- 



(iii) 

K 



E, 



(ii) 



(17) 
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= 2h (cos K x + cos K y ) + r (cos 3K X + cos 3K y ) , 

the results being summarized in Table ^ at the fourth 
order of a t/U expansion with: 



E, 



(IV) 



t 2 r 4 
0.9023C/ - 208.9— - 682883- 



(18) 



U U 3 
and h = 1000t 3 /U 2 , r = 3320t 3 /C/ 2 and s = 61926i 4 /[/ 3 . 



K- 6/2vr 




(0,0) 


E { 1V> - Ah + 2r + As 


(1,0) 


E ( lv) -Sh + s 


(1.1) 


E { 1V) -2h-2r-2s 


(2,0) 


E[' V) -h + 2r + s 


(2,1) 


El' v) -2s 


(2,2) 


E[' V) +2h + 2r~2s 


(3,0) 


E^ v >+4s 


(3,1) 


E^ v) +h-2r + s 


(3,2) 


E { IV) +3h + s 


(3,3) 


E[' V) +4h-2r + As 



Table 1. The L 2 lowest energies given by a t/U expansion up 
to order four. The other K-states follow from the x-y, and the 
K — > K symmetries. 



One can see in Fig. that this t/U expansion gives 
the exact momenta K of the 36 first states above a last 
level crossing at r s = 196. For L = 6, one has no GS 
level crossing, having the same four total momenta K = 
27r/6(±l,±l) for the four GSs when t — > and when 
U -> 0. 

This absence of GS level crossing is a property of the 
case TV = L/2. If L/2 > N = 3, the GS momentum is 
K = when t — > 0, and a GS level crossing takes place 



as r s increases. This correlated lattice regime in the limit 
L — > oo is discussed in App. [5] 

Let us now consider one of the GS wave functions of 
momentum K, for instance K = 27r/6(1,1). When t — > 0, 
it reads 

|fo(K)) = l^ e ^£)j|0), 



(19) 



and corrections are given at the leading order (first order) 
of the t/U expansion by: 

a E n — E n 



(20) 

Only twelve t — > eigenstates |!^ Q '(K)) are directly cou- 
pled to |tfo(K)) by a single hop (coupling term — t). This 
gives: 



12 



^ (/) (K)) = |^(K))+£ 

a' = l 

(0) 



t 



AE n , 



hMK)), (21) 



where AE a > — E a i — E n . The t/U expansion of the av- 
erage (!?o(K,r s ) |/|!?o(K, r s )) of an observable / which 
takes definite values over the states |<? a /(K)) are given at 
leading order (a (t/U) 2 ) by: 



</} (//) = 
where 



(ii) 12 / . 

q' = 1 



\AE a , 



f(*a'l (22) 



a" 



AE a . 



(23) 



(77) 12 
a'=l 

This can be used to calculate the n t/l moment of the 
different inter-particle spacings in the correlated lattice 
limit. Each state 



c\c\c\ |0) 



(24) 




50 100 150 200 250 300 



0.01 



0.001 




1000 



Fig. 5. The 36 lowest energies as a function of r s for all possible Fig. 6. The relative fluctuations u r of the three inter-particle 
total momenta K. Above the last level crossing indicated by spacings r m in (solid line) r; n t (dashed-dotted line) and r max 
the arrow the momenta of the 36 first states are given by the (dashed line) as a function of r s . The thin dotted lines give the 



tj [/-expansion (Table^. 



perturbative t/U-dec&ys given after Eg. 1261 
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is characterized by three inter-particle spacings 
r m in(ijl) < r int (ijl) < r max (ijl). Taking for / 



f(r n ) = J2m) r n W) (Mil 

ijl 



(25) 



one can calculate the averages (r) , the variances ( r 2 ) — (7-) " 
and the relative fluctuations 



u r 



<{r 2 



(26) 



of Train-, Tint and r max in the correlated lattice limit. 

The three relative fluctuations u r are shown in Fig. |SJ 
as a function of r s . One can see that above r*, they co- 
incide with the behaviors u r (r m in) ~ 5A9t/U, u r (ri n t) ~ 
AAlt/U and u r (r max ) « 12.34t/£/ given by the leading 
order of the t/U expansion. 



0.1 - 




_i 1 I 1 1 I 1 1 I 1 1 



0.01 



0.1 



10 



100 



Fig. 8. The relative fluctuations u r of the three distances r m i„, 
Tint and r max as a function of r s . The dotted lines give the HF 
values, which coincide to the actual behaviors up to rf F ~ 1. 



6 Breakdowns of the perturbative expansions 

When a perturbation is weak, the eigenstates remain local- 
ized in the vicinity of the unperturbed states. Increasing 
the perturbation delocalizes the eigenstates in the unper- 
turbed eigenbasis, yielding a crossover fTI] from a weak 
perturbative mixing of the unperturbed eigenstates (Rabi 
oscillations) towards an effective golden-rule decay. Above 
the derealization threshold, an expansion around the un- 
perturbed eigenbasis does not make sense. This dereal- 
ization threshold is given by a general criterion discussed 
in different contexts: onset of quantum chaos in a many 
body spectrum |11I12I13| . quasi-particle lifetime and de- 
localization in Fock space |14ll5| h interaction induced 



T 1 1 1 — I I I I J 1 1 1 1 — I I I I J 1 1 1 1 — I I I I 




1 I 1 . / I 1 I 1 

1 10 100 1000 



Fig. 7. First energy spacing AE(r a ) (dotted line) with its two 
limits (solid lines - 2t for U -» and 0.025(7 for t -> 0), and 
spreading widths Fp^s) (dashed line) and IV (r s ) (dashed- 
dotted line) calculated from the local densities of states Qf(E) 
and Qw(E) respectively. L — 6 and K = (2ir/6, 2ir/6). The 
arrows indicate the limits rf and r* between which the per- 
turbative expansions are inappropriate. 



Q_ 




0.1 1 10 100 1000 



Fig. 9. The GS projection | (>f r o(r s )|<f r o(r' i , = 0)} j 2 over the non- 
interacting GS (or HF-GS) as a function of r s . Notice the 
breakdown of the HF-behavior at rf F ~ 1. 

thermalization |16) : Quantum ergodicity occurs when the 
perturbation matrix element (i\P\f) between an unper- 
turbed eigenstate \i) to the "first generation" of unper- 
turbed eigenstates |/) directly coupled to it by the per- 
turbation is of the order of their level spacing Ef — Ef. 

(i\P\f)^Ei-E f . (27) 

Using this general criterion for the three particle GS, one 
can estimate the range of intermediate ratios r s where the 
GS cannot be simply described neither in terms of the low 
energy weak coupling eigenstates, nor in terms of the low 
energy strong coupling eigenstates. 

6.1 Limit of the weak coupling [//£-expansion 

When U/t — * 0, one can use the Hartree-Fock (HF) 
mean-field approximation for the ground state. Since the 



Z. A. Nemeth and J.-L. Pichard: Andreev-Lifshitz supersolid II 



7 



1 r 



0.1 r 



0.01 - * 



0.001 



|V(0,0) 
|V(27t/L,0) 
|V(27t/L,27t/L) 




_i i i i_ 



_i i i i i i i i_ 



0.05 
0.04 
0.03 
0.02 
0.01 




-7t / 


1 II 1 ,1 1 

-4t; -t 


ii i in ii | ■ A 

\ 21;' 5t/\ 


j 


1 
1 


! i f \ 

' i 1 -i 

' 1 




Tl 




I r 


3 ' ' i 






//5t ; 




/ \ 












/ -8t/ 




/ r'7 ~\'\- 















10 



100 



10 



100 



Fig. 10. The dimensionless interaction matrix elements | V(q)| 
(see Eq.QJ as a function of L for different values of q, exhibiting 
the asymptotic 1/L decay. 

charge distribution is uniform without disorder and with 
periodic BCs, this approximation becomes trivial. Start- 
ing from plane wave states of uniform density, the HF- 
approximation consists in studying again a single particle 
in a uniform potential. Therefore, the GSs for U — re- 
main the self-consistent states of the HF-approximation, 
while the HF-eigenenergies are given by the first order 
pertubative expressions (see Eas. HTl and IT2*|) . The range 
of validity of the HF-approximation can be seen in Fig. 
El Fig.|Hland Fig. for various observables. As one varies 
r s and for L = 6, Fig. |2] gives the GS-energies, for the 
total momenta K = 2%/L(l, 1) and K = 27r/L(0, 0) re- 
spectively. The relative fluctuations u r (U) of the three 
inter-particle spacings and the projection 

2 



Phf(U) = \{%(U = 0)|!Pb(E0>r 



(28) 



of the actual GSs onto the HF-GSs are given in Fig.|Hland 
Fig.Elrespectively, for a total momentum K = 27r/6(l, 1). 
These figures show that the HF-approximation remains 
valid up to a value r^ F w 1. 

This numerical value for rf F is close to the value 
r F yielded by the general criterion (Eq. I27fl giving the 
crossover from a weak perturbative mixing of the unper- 
turbed states towards an effective golden-rule decay (de- 
localization threshold in Fock space). 

Let us first give a qualitative estimate valid for ar- 
bitrary values of L and N, where the values of U, t, N 
and L appear through the expected dimensionless ratio 
r s . Assuming L sufficiently large for having cos(27r/L) w 
1 — 2{tt/L) 2 , the energy spacing AE = E\ — E between 
the GS and the first excitations for U = reads: 



AE = 2-Kt ^2 



kl - k%) » 2vrf (2k F q) 



WN 



(29) 



since kp ~ vTV /L and q ~ 1/L. The interaction matrix 
element directly coupling these states reads: 

U 
L' 



\H int (<i)\l) 



(30) 



Fig. 11. GS projection | (S>b(r 8 )|!P„(r 8 = 0)) | 2 onto the 
non-interacting eigenspaces of increasing energies E a — 
—St,—7t,... and total momentum K = (2tt/L, 2tt/L) as a 
function of r s . Notice the behaviors of the projections onto 
the \\P a (r a = 0)} which do not contribute to |^o(r s — * oo)) 
(E a = —8t, 



as shown in Fig. 111)1 for a large enough L. The criterion 
(Eq.EZ) AE « (0\H mt {q)\l) gives: 



const. 



const. 



(31) 



For an exact determination of rf on 6 x 6 lattice, 
we have numerically studied the local density of states 
(LDOS) H8, 

q f {E) = Y,\{*o(U = 0) I #i(U))\ 2 6(E -Ei(U)) (32) 



of the non-interacting GS in the eigenbasis with interac- 
tion. This LDOS is a distribution of width _Tp, where h/Tp 
gives the lifetime of the U — states when U is turned 
on. Using Fermi golden-rule for estimating fp: 



r F = 2tt ^2 1 (0|ff<„t|l) \ 2 9(Ei) - U 2 \V(a)f 



1 



El — Eq 1 

(33) 

and since the density of states g(E\) » 1/AE, criterion 
(l2*7jl corresponds to 



r F w AE. 



(34) 



For K = 27r/6,27r/6), the GS of energy -lOt is directly 
coupled to two states of energy — 8t by a matrix element 
2[/(V r (27r/6,27r/6) - V(2tt/6,0)) = 0.091J7, and the con- 
dition 

r F =4J 0m ™r*AE = 2t, (35) 

is satisfied when rf ~ 6. We have also calculated r F di- 
rectly from q f (E), and one obtains again r = AE at 
rf w 6, as shown in Fig. One can notice that 
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The ratio rf (breakdown of a U/t expansion) is larger 
than the ratio r^ F fa 1, where the HF-behaviors cease 
to be valid (zero order approximation for the wave 
function, first order approximation for the energies). 
Fig. ^2 gives the GS projection over the non- 
interacting eigenspaces of increasing energies. The pro- 
jection onto the first excitations |1) of energy — 8t is 
particularly interesting. These states are directly cou- 
pled by the interaction to the GS for U = 0, but are or- 
thogonal to the GS when U — > oo. One can see that the 
increase of this projection driven by the U/t expansion 
ceases precisely at rf. This is the first manifestation 
of the correlated limit as U increases. 



6.2 Limit of the t/U correlated lattice expansion 

We use the criterion (Eq. 127(1 for determining the value r* 
under which the correlated lattice expansion breaks down. 



The matrix element 



\Hkin\l) 



t is equal to the first 



corresponding energy spacing AE for U = U* 
satisfies 

t = AE » \/2— + 0(£~ 3 ), 



L 2 



where U* 



(36) 



which yields the threshold: 



U* 

Tl 2 



const 



const ■ L 



(37) 



When r s < r* oc L 3 , the GS becomes delocalized in the 
r s — > oo eigenbasis and the t/U expansion breaks down. 
This gives the same L-dependence for r* than in Sec. |5J 
From the local density of states 

qw{E)=Y^\{MU = oo) \MU))\ 2 5{E-E i ), (38) 



of the GS \\Pq(U = oo)) when U — > oo in the eigenbasis 
with interaction, we have calculated the spreading width 
ryy of \&o(U = oo)) when one turns on the kinetic hopping 
term t. One can see in Fig. [7| that /V = AE for r* fa 
180. Above r*, the three inter-particle GS spacings are 
correctly given by the {/-expansion (see Fig. |SJ) and the 
ordering of the L 2 low energy levels correspond to those 
given in Table ^ the last level crossing being shown in 
Fig. at r s fa r*. 



7 Partially melted Wigner molecule near the 
correlated lattice limit 



a square with a 'vacancy' at one of its corners. There- 
fore, beside the rigid translation of the triangle (hop- 
ping term h oc t 3 /U 2 ), there is another possible third 
order process: the tunneling of the vacancy (hopping 
term r oc t 3 /U 2 ). As r s decreases, the GS may have ad- 
vantage to delocalize this vacancy to reduce its kinetic 
energy. This will be the "Andreev-Lifshitz supersolid" 
in an ultimate mesoscopic limit. 

If one looks at Fig. |H| we see that the smallest inter- 
particle spacing r m i n begins to fluctuate according to 
the lattice i/{/-expansion above r s w 40 while the 
largest one r max enters into this lattice regime at a 
higher value r s fa r* fa 180. The behaviors of r m i n 
and r max suggest us that one has for 40 < r s < 180 
a partially melted Wigner molecule (PMWM). This 
PMWM is composed of a rigid two particle Wigner 
molecule (2PWM) with oscillations around the equilib- 
rium positions smaller than the lattice spacing a (con- 



sistent with the correlated lattice behavior of 



0. 



while the third particle remains more delocalized, ex- 
plaining the absence of a correlated lattice behavior of 

'max • 

One can see the derealization of the third particle in 
the direction parallel to the 2PWM formed by the two 



O--O--G--O--0--O 

K-0 





We begin to study the GS in the non-perturbative regime 
rf < r s < r*. We first introduce a simple ansatz proposed 
in Ref. [7] and which turns out to describe the GS for 
r s fa 40. The idea of this ansatz can be given by three 
observations: 

— As shown in Fig. [3] and Fig. 0] when r s is large, the 
three particles form a triangle, which can be seen as 



Fig. 12. Above: Scheme of a a:-oriented PMWM: one particle 
is totally delocalized in a;-direction with a momentum k, but 
remains localized in the {/-direction at a distance L/2 — 3 
from the two other particles which form a 2PWM. Below: 
The GS projection amplitude ip(x,y) (see Eg. 1391 at r s = 40. 
ip(x,y) is not totally uniform in the a:-direction, since there is 
the contribution of the y-oriented symmetric PMWM. 
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Fig. 13. The GS-projections for K = 27r/6(l, 1) over different 
x-oriented PMWMs ^(K): solid (« = 0), dash-dotted (k = 
4n/6) dotted (k = 87r/6) lines respectively. The upper thin 
dashed line gives the total GS-projection over x «-» y symmetric 
l#«=o(K)> 



Fig. 14. GS-projection over the PMWM-ansatz # K=0 (K) 
when the oscillatory motions of the PMWMs are included up 
to 0th (dashed), 1st (dashed-dotted) and 2nd (dotted) orders 
of a lattice t/!7-expansion. 



others, in Fig.^J where the GS projection amplitude 



ip(x,y) = (ffb(r») 
is given for r s =40. 



C (l,l) C (4,l) C (x,») 



(39) — 



Let us consider a x-oriented state where two particles 
are fixed, while the third particle is delocalized along one 
line with a momentum shown in Fig. 1121 



C j C j+(3,0) C lj„+3 



|0> . 



(40) 



To satisfy translational symmetry, one defines the x- 
oriented PMWM of total momentum K by 



i(Kj— Kj x 



N 



c ] c ] + (3,o)4 JB+ 3)l°)- ( 41 ) 



0.1 




10 



100 



Fig. 15. Root mean square Ar of the three inter-particle spac- 
ings rmi n , Virvt and r max as a function of r s : actual behaviors 
(thick lines) and behaviors given by the 2nd order PMWM- 
ansatz (dotted lines) . The arrows underline the close agreement 
for r a « 40. 



Then, we combine the x-oriented PMWM with its sym- 
metric y-oriented counterpart: 



l^(K)> 



^(K)) 



(42) 



for satisfying the x <-* y symmetry. Moreover, k can only 
take three quantized values (0, 47r/6, 87r/6) for L = 6. 

One can follow the delocalization of the third par- 
ticle in Fig. El When r s — > oo, the GS-projections 
|(^(K) \%{r s ))\ 2 over the x-oriented PMWM states of 
K = 0, 4-7I-/6, 8n/6 become equal. This corresponds to a 
total localization of the third particle on the line parallel 
to the 2PWM, such that the three particles form rigid tri- 
angles. When r s < 40, only the contribution with k = 
is non-zero: the third particle becomes fully delocalized. 



Moreover, if we take the x <-> y symmetric combination, 
we can see that we describe more than 90 % of the the 
real GS when r s w 100 < r*. If r s is further decreases, 
our ansatz have to include some oscillatory motions of 
the 2PWMs around equilibrium. If these oscillations are 
smaller than the lattice spacing, it is enough to use a 
similar lattice t/£/-expansion as in l|2()|l . for the 2PWMs 
only. Fig. 1141 shows the improvement of the GS-projection 
\ 2 

<f K (K) | &o(U)) when the PMWM oscillations are in- 
cluded. The GS-projection over the expanded ansatz ex- 
ceeds now 95 % at a smaller value of r s ps 40 if the t/U- 
expansion is extended up to the 2nd order. 

Fig. 1151 shows how the PMWM-ansatz, with 2nd or- 
der oscillatory motions of the 2PWMs, describes the root 
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mean square of the distribution of the three inter-particle 
spacings r rnin , r int and r max around r s ps 40. 

Above r* w 180, the "triangles" are formed. Below 
r*, the gradual delocalization of one particle parallel to 
a remaining 2PWM with increasing oscillations begins. 
Around r s ss 40, this delocalization is completed with a 
momentum k = 0. This gives a mixed state, having both 
the properties of a "solid" and of a "liquid" . However this 
mixed state exhibits lattice effects and cannot be the final 
step of the melting process. 



8 Partially melted Wigner molecule near the 
Fermi limit 

When the interaction is further decreased, two things have 
to be taken into account. 

— One concerns the "solid" : the amplitudes of the oscil- 
latory motions begin to exceed the lattice spacing, and 
cannot be described by the lattice t/U expansion. 

— The other concerns the "liquid" : the Coulomb repul- 
sion cannot be strong enough to maintain the delocal- 
izcd particle parallel to the correlated pair. 

To describe the next step of melting, in the range 
rf < r s < 40, let us consider the Fermi limit and the 
picture |S] suggested by Bouchaud et al of a system of un- 
paired particles with a reduced Fermi energy co-existing 
with strongly paired, nearly solid assembly. For the ul- 
timate limit N — 3, this simply means a single corre- 
lated pair co-existing with a third particle remaining in its 
non-interacting state. The third particle roughly remains 
in a plane wave state of uniform density (momentum 
k 3 = (0,0), 271-/6(1,0) or 2tt/6(0, 1) if K = 2tt/6(1,1)), 
and provides a uniform background for the correlated pair 
formed by the two other particles. The occurrence of a 
correlated pair is suggested by the breakdown of the U/t- 
expansion above rf on one side, by the previous PMWM 
ansatz on the other side. The conservation of the total 
momentum allows us to order these 2+1 particle systems. 
Different combinations, given in Tabled are possible. 



0.01 



_i i ' i i ' i i ' 



0.1 



10 

u 
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Fig. 16. Above: N = 2 particles: Relative fluctuations u r (U) 
of the pair GS (solid line - asymptotic spacing a = (3,3)) 
and first excited state (dashed line - a = (3,2) or (2,3)). The 
dotted lines give the three regimes for the pair: (Fermi liquid 
- continuous Wigner solid - correlated lattice solid). Below: 
N = 3 particles: GS-projection Po{U) (dashed line) on the 
non-interacting GS, and GS-projection (Pa(U), solid line) and 
(P^ ! \U), dotted line) onto the four states given in table 
the pair wave function being calculated exactly or using the 
2nd order i/C/-expansion respectively. 





k3 Q 


K 2a 






2tt/6(0, 1) 


2tt/6(1,0) 


(3,3) 




2tt/6(1,0) 


2tt/6(0, 1) 


(3,3) 


*3 


2tt/6(0,0) 


2jt/6(1, 1) 


(3,2) 


$4 


2tt/6(0,0) 


2tt/6(1,1) 


(2,3) 



Table 2. Four different combinations of an unpaired particle 
of momentum k.s a co-existing with a correlated pair of mo- 
mentum K.2a- + k3 Q = K = 27r/6(l, 1) and a„ is the pair 
inter-particle spacing when U — ► oo. 



In this table, a Q is the asymptotic value of the inter- 
particle spacing of the pair, if we follow the level to the 
limit U — ► oo, the asymptotic pair wave function becoming 
|0). However, for K2 = 27r/6(l, 1), this pair wave 



function with a = (3,3) is zero. In this case, the pair 
GS has a twofold degeneracy, with a = (3, 2) and a = 
(2,3). This is why we have four combinations instead of 
three, when we keep one particle out of three in its non- 
interacting wave function. 

Let us begin to consider those four states \& a )> assum- 
ing that the pair behaves as rigidly as in the limit r s — > 00: 



10) 



(43) 



Let us first allow oscillatory motion of the pair using the 
2nd order lattice £/?7-expansion. In Fig. El one can see 
the total GS-projection 



(44) 
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onto these four states. As U decreases, P^^ {U) increases, 
but saturates below U « 20. This saturation tells us that 
the lattice t/U perturbation theory is no longer a suit- 
able tool for describing the pair. To improve the GS de- 
scription, we numerically calculate the exact wave func- 
tion \$K. 2a (JJ)) °f the pairs: 



UP 



K, 



,(P))= E 



^(^4410); (45) 



k 1 +k 2 =K 2Q 



to have the following ansatz wave function for each com- 
binations: 



\ 



E ^(U)clcl 



\ k l. k 2 
\k 1 +k 2 =K 2Q 



cLjO)- (46) 



The total GS-projection 



Pa(U) = 52\(¥«(U)\MU))\ 

(a) 



(47) 



over the subspace spanned by the four ansatz states is 
given by a solid line in Fig. 1161 fbelowh For comparison, 
we have also plotted the GS-projection 



p (u) = \(Mu = o) \MU))\ 



(48) 



onto the the non-interacting GS (dashed line) of same K. 
If one uses unpaired fermions co-existing with correlated 
pairs instead of the non-interacting GS, one substantially 
improves the description of the intermediate GS above 
rf. For U « r s = 10, the projection Pa{U) exceeds 90 
%. In Fig. El on e can see for the same value of U the 
relative fluctuation u r of the corresponding two particle 
system, both for the a a — (3, 3) and (3, 2)-states. This 
study of the case N — 2 shows that the pairs are nei- 
ther in their solid lattice regime, nor in their liquid non- 
interacting limit when U ~ 10, but in their continuous 
Wigner regime. 



9 Andreev-Lifshitz supersolid for intermediate 
couplings? 

With this second ansatz, our understanding of the differ- 
ent steps of the melting of the three particle system on 
a 6 x 6 lattice is achieved. Above r* w 180, we have de- 
localized correlated triangles and a lattice t/U expansion 
is sufficient for describing the L 2 low energy states. For 
the GS, the melting begins with the partial derealization 
of one particle in the x or y directions down to r s « 40. 
This process is the natural extension of the vacancy tun- 
neling which is the "softest" degree of freedom of the cor- 
related triangle, if one views a triangle as a square with 
a vacancy. This derealization of a single particle is first 
one dimensional, parallel to the remaining 2PWMs. For 

10 < r s < 40 there is a crossover in the GS-structure. One 
particle becomes now totally free while the other two form 



a floppy pair with oscillatory motions. One can view this 
supersolid GS as a natural extension of the mean field HF 
state. Instead of having one particle in the mean field of 
the others, one has a pair in the field of the third particle. 
This is reminiscent of the BCS ansatz with a fixed num- 



ber of fermions proposed in Ref. Below r s 



, the 



remaining pair melts and one recovers the HF mean field 
limit. 



As explained in Sec.|3|for N = 2, above the first thresh- 
old rf , one has a continuous Wigner regime (oscillations 
of the molecule around the equilibrium positions exceed- 
ing the lattice spacing), before having important lattice 
effects at a second threshold r* cx L 3 . As clear in Fig. 
^ and Fig. L = 6 is too small to have the expected 
correlated three particle Wigner molecule in a continu- 
ous regime, without lattice effects. L — 6 is just large 
enough for having a partially correlated supersolid regime 
free of significant lattice effects. However, the exact di- 
agonalization study of N = 3 spinless fermions can be 
extended to larger values of L. Such scaling analysis is 
in progress. The first results confirm the general picture 
emerging from this detailed study of the L = 6 system, 
and support, at least for an ultimate mesoscopic limit, 
the possibility proposed by Andreev and Lifshitz for the 
thermodynamic limit: a quantum crystal may have delo- 
calized defects without melting, the number of sites of the 
crystalline array being smaller than the total number of 
particles. 



On one side, one cannot exclude, as already pointed 
out in Ref. fj], that a supersolid regime is favored by the 
chosen geometry, because of the number of particles and 
underlying square lattice, but not favored at all in the 
continuous limit, which has not a square symmetry but a 
spontaneously broken hexagonal symmetry. On the other 
side, such a supersolid regime is not considered in the 
quantum Monte-Carlo (QMC) studies, since these meth- 
ods rely on certain assumptions on the GS nodal struc- 
tures. We study in App. El those nodal structures in our 
small lattice model, which turn out to be very complex for 
intermediate values of r s , and cannot be approximated by 
the nodal structures of the weak or the strong coupling 
limits. In App. [Dl the GS occupation numbers in the re- 
ciprocal lattice are given when r s increases, and exhibit 
a hybrid structure in the mesoscopic supersolid regimes, 
the usual spreading of the "Fermi sea" being accompanied 
by the gradual emergence of the Fourier spectrum of the 
"Wigner solid" . 
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A Harmonic oscillatory motion of a two 
particle molecule 

In this appendix, the oscillatory motion of a two particle 
Wigner molecule on an empty L x L lattice with periodic 
BCs is studied when rf < r s < r*. This corresponds to 
a Coulomb repulsion which is strong enough for forming 
a Wigner molecule, but not strong enough for restricting 
the oscillatory motion of the particles around equilibrium 
to scales of order of the lattice spacing a (correlated lat- 
tice regime). As earlier noticed, the distance between two 
sites on a periodic lattice defined in Eq.|3] yields a cusp of 
the Coulomb repulsion at the equilibrium positions of the 
Wigner molecule. The role of this cusp is negligible in the 
Fermi limit, but not in the Wigner limit, the oscillatory 
motion of the particles around equilibrium being located 
at the cusp. To avoid this complication, we smear the long 
range part of the Coulomb repulsion, taking for the pair- 
wise repulsion between two particles separated by r in the 
continuum limit of a system of size D 



V c (r) 



e 2 7r 



D 



(49) 




Fig. 17. Re-definition of the domain of variation for R and 
r in one dimension. The gray area is the original domain < 
xi , X2 < D with periodic BCs. Moving the part where R < 
to the upper triangles, the new domain becomes < R, r < D 
with periodic BCs. The bottom left empty circle, which is the 
x\ <-> X2 exchange counterpart of the solid circle, goes to the 
empty circle in the upper right triangle. 



D 



instead of the previous repulsion V(r) = e 2 /|r| with |r| 
defined by Eq.[3]for a square with periodic BCs. The cor- 
responding two particle Hamiltonian 



h c = -—{v\ + vD + v c {v) 

2m 

is the sum of two decoupled terms 

H c = H C (R) + H c (r) 



(50) 



(51) 



when one uses the center of mass R = (r^ +r 2 )/2 and rel- 
ative separation r = ri — r 2 coordinates. The Hamiltonian 
for the center of mass motion is given by 



h c (h) 



4m 



R 



while the Hamiltonian for the relative motion reads 



H c {r) = V 2 + t/(r). 

m 



(52) 



(53) 



It is convenient to redefine the domain where R and r 
vary such that it becomes again a 2c? torus. Fig.ll7lshows 
how it can be done in one dimension. This change of do- 
main modifies the usual symmetry requirement 

&(R,r) = y{x x ,x 2 ) = -!?(x 2 ,a?i) = -W(R,-r), (54) 

for spinless fermions, since ^j(R, — r) — > &(R + D/2, — r + 
D). In two dimensions, the symmetry requirement for 
spinless fermions after this change of domain becomes: 

#(R, r) = -*{r + {E,EY-t+ (D, £>)) ■ (55) 



The eigenstates of H c = H C (R) + H c (r) take the form 

1 



!P(R,r) = -e u "ty(r) 



(56) 



the symmetry requirement for ip( r ) being: 

W(A D)-r), for §(K X + K y ) even; 



V>(r) = 



— r 



for §(K X + K y ) odd. 



(57) 

When a sufficient Coulomb repulsion yields small os- 
cillation of the inter-particle spacing around its largest 
possible value, the Coulomb repulsion can be expanded 
and one gets for the relative motion of the two particles a 
2d harmonic oscillator Hamiltonian: 



HJr 



e 2 7r 3 



DV2 D 3 V32 
For a 2d harmonic oscillator 

h 2 









t) 



2m 1 



(58) 



(59) 



the symmetric and antisymmetric GSs of energies Eqs = 
%u> and Eqa — 2hu are given by: 



1 



exp 



and 



,T, I \ T ^ 

^oA(r) = _ exp 



3/2 

o y 



2/ 2 
zt (i 



2/ 2 



(60) 



(61) 
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Fig. 18. GS energy Eo for two particles in a 12 x 12 lattice 
as a function of U. We have plotted for convenience E = Eq — 
0.185(7 — 0.137 (taking i = 1) for checking the contribution 
~ (7 1//2 at intermediate U of Eq. 1671 The analytical estimate 
(0.15907, dotted line) is valid for 1 < U < 1000, defining 
the continuous Wigner regime where the relative fluctuation 
u r oc r7 1//4 (Fig0. 



respectively. The length Iq = y/h/ (rruJ) becomes in our 
case 



h 



(62) 



Discretizing the continuous Hamiltonian on L x L lat- 
tice, one gets 



Hl = -t^ c\ C y + Wt 
<JJ'> 



Utt 



L\ /sin 2 (J-* 



sm 



2 (jyl-jy2)n 

D 



(63) 



where L = D/a, U = e 2 /a and t = % /2ma 2 . The constant 
term ANt comes from the discretization of the Laplacians. 
For Hi, the characteristic length Iq becomes 



(64) 
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Fig. 19. N = n — 1 spinless fermions on a square lattice of 
size L = 3n for n — 4. One of the L 2 square Wigner lattices 
with a vacancy of energy _Eq ' (solid circles) and one virtual 
state of energy E[ (empty circle) contributing to the vacancy 
dynamics are indicated. 



where 



TlLO = 



Uh 2 l e 2 7r 3 



(66) 



This energy becomes for the corresponding lattice Hamil- 
tonian Hf. 



En 



2tt 2 Utt 

-TT t + 7= 

L 2 LV2 




(67) 



A numerical check of this expression using a 12 x 12 lattice 
model is shown in Fig. ED Dividing Eq by U 2 /t, one gets 



Etit 
U 2 



4r? 



r 3/4 



4r s 



4r. 



3/2" 



(68) 



This expansion is similar to the original expression given 
by Wigner |!5] for the strong coupling limit, the energy 
being measured in Rydberg units. The first term gives the 
kinetic energy of the center of mass (~ r~ 2 ), the second 
is the electrostatic energy at equilibrium (~ rj 1 ) and the 
third term comes from the oscillations of the inter-particle 
spacing around equilibrium ( 



-3/2, 



Since lo oc r s 1 ' i D in our model, one gets (r) = \/2D 

for the average inter-particle spacing and Ar oc r s 1 ^ 4 D 
for the width of its distribution. Therefore, the ratio u r of 

— 1/4 

these two quantities decays as r s ' . This corresponds to 
the behaviors shown in Fig. ^ when rf < r s < r*. 

For a GS of total momentum K = 27r/D(l,0), the 
condition (|57|l yields the symmetric GS for H c (r). For the 
two particle GS of H C (H) + H c (r), one gets an energy 

h 2 e 2 n 

E a = — K 2 + L + hu, (65) 

4m DV2 



B Correlated lattice limit when L — > oo 

We study two limits which can be easily described by the 
t/U lattice perturbation theory when t —> and L — > oo. 

We keep N = 3 spinless fermions in the first limit. In 
this case, the hopping term ~ h characterizing the rigid 
translation of the molecule remains of order t 3 /U 2 , while 
the hopping term ~ r characterizing a single particle hop 
over L/2 sites and coupling triangles of different orienta- 
tions becomes of order t(t/i7) i / 2_1 . When L — > oo, only 
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the rigid translation of the triangle matters, and the ef- 
fective Hamiltonian fEa. I16j) reads: 



H 



(in) 

eff 



E^« 







<j,j'> 



(69) 



where -ft (Dj) creates (annihilates) a triangle defined by 
Eq. El replacing 3 by L/2. The resulting energies are: 



E 



(in) 

K 



E, 



in) 



2h (cos K x + cos K y ) , 



(70) 



yielding a GS total momentum K — (0, 0). 

We take N = n 2 — l spinless fermions and a size L = 3n 
in the second limit. N — n 2 gives a uniform filling factor 
v = 1 /9 and a square Wigner lattice which is commensu- 
rate with the assumed lattice. Taking a single particle out 
n 2 will create a single vacancy in the square Wigner lat- 
tice, as shown in Fig. ED In this second limit, this is now 
the rigid translation of the Wigner lattice which becomes 
negligible in the thermodynamic limit, while the hopping 
term r characterizing the propagation of the vacancy re- 
mains oc t 3 /U 2 . The effective Hamiltonian fEa. I16fl takes 
the form: 



H 



(in) 

eff 



E 

(J J') 3 



D\Dy. 



(71) 



There is only a set of virtual states of energy E[ 0) which 
contribute at the lowest order to the propagation of the 
vacancy (see Fig. 119(1 . making simple to calculate r. After 
Fourier transformation, the eigenenergies of the effective 
Hamiltonian are given by: 



E 



(in) 

K 



= E, 



(ii) 



2t 3 (cos 3K X + cos 3K V 



(4 



(0) 



p(0) 



(72) 



and corresponds to the spectrum of a single particle on 
the assumed square lattice with third nearest neighbor 



hopping. Of course, this one vacancy dynamics does not 
totally remove the L 2 GS degeneracy of the limit t — > 0, 
but makes very simple the study of charge propagation in 
this highly correlated many particle system. 



C Nodal structure of the three particle 
system 

The quantum Monte-Carlo (QMC) methods are the most 
powerful tools for studying large many-body systems |17j . 
However, the study of the ground state of fcrmionic sys- 
tems suffers from the well known "sign problem". One 
way to avoid the negative weights that would be other- 
wise generated by antisymmetric states is the fixed node 
approximation |17| . The fixed node GS energy is then an 
upper bound to the exact GS energy. The nodal structure 
of the liquid limit is given by a simple Slater determinant 
of plane waves, and of localized orbitals for the solid limit. 
To know the exact nodal structure for intermediate cou- 
plings is not an easy problem 18 . In this appendix, we 
study the nodal structure of three spinless fermions on a 
6x6 lattice. 

Previously, we have considered K-eigenstates, for an 
interacting system which is invariant under lattice trans- 
lations. The Hamiltonian being invariant under time- 
reversal symmetry, one first define eigenvectors with real 
components in the site basis. For this purpose, we combine 
the K-eigenvector 

|*k(r,)>= E *>«to*»( r >)<i 1 <L<L 10), (73) 

k l . k 2 k 3 
k 1 +k 2 +k 3 =K 

with its time reversed conjugate —K-eigenvector: 



,k2 ,1^3 
k 1 +k 2 +k 3 = K 



(74) 



900 I 




0.1 1 10 100 1000 



Fig. 20. The number of negative components -/V„(0) as a 
function of r 3 . The solid line is obtained with a threshold £ = 
1CF 12 . The dotted line is a guide to the eyes. 



1200 i ■ ■ i i 




0.1 1 10 100 1000 



Fig. 21. The number of negative components N n (oo) as a 
function of r 3 . The solid line is obtained with a threshold £ = 
10~ lj . The dotted line is a guide to the eyes. 
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possible values 1/9 (filled circles) and (empty circles). 




Fig. 22. Occupation numbers Pk on the reciprocal lattice 
(2nk x /L, 2-irky/L) for r s = 0.1 (upper figure), r s — 10 (middle 
figure) and r s — 40 (lower figure) . 

Since the tf'ki,k 2 .k 3 {r s ) are real, the combination 

*h,h,h( r s) = -j= (0\c h c h c h [\^ K (r s )) + |#-_ K (r,)}) 

(75) 

is indeed real in the site basis. 



We want to compare the nodal structure of this real 
GS at intermediate values of r s to the two limiting nodal 
structures, characterizing the GS either in the limit r s — ► 
or r, — > oo. Wc define the components of a vector 
X(r a ,lim) by 

Xh, h,h( r s,Hm) =%j 3 j 3 (r a ) •%j 2 j s (iim); (76) 

where ^(lim) is the corresponding limiting GS, and we 
count the number N n (lim) of negative components of \. 
If the nodal structures of the limit and at r s are identical, 
N n (lim) = 0. However, when Xji ,h,h ( r s> Hm) is almost 
zero, its sign becomes undefined due to numerical preci- 
sion. This is why we ignore all the components of x below a 
given threshold £. The behavior of N n (0) (N n (oo)) where 
&{lim) is a liquid GS W{r s = 0.1) (solid GS &{r a = 10 3 )) 
is shown in Fig. EDI (Fig. EJ. 

As one can see, the nodal structure does not exhibit a 
sharp transition between the two limits, but a crossover 
with complex intermediate behaviors. Notably, there are 
some plateau values suggesting some constant nodal struc- 
ture around the intermediate values of r s where we observe 
the PMWM states. This illustrates the difficulty to use a 
fixed node Monte-Carlo method for describing the inter- 
mediate GS on a lattice. 



16 



Z. A. Nemeth and J.-L. Pichard: Andreev-Lifshitz supersolid II 



D Occupation numbers in k-space 

It is not only interesting to know how the system oc- 
cupies the real lattice, but also the reciprocal lattice. 
From the GS wave function ^(js) of total momentum 
K = 27r/6(l,l), we have calculated the occupation num- 
bers 

flc(r.) = (*b(r.)|4dk|*b(r.)), (77) 

in k-space. Fig. I22land Fig. 1231 give plots of those number 
in the reciprocal lattice for increasing values of r s . When 
r s — > 0, only three k-states are occupied (Tig. 12211 . When 
r s — ► oo, one gets a simple pattern (Fig. I23fl for K = 
27r/6(l, 1), which can be analytically obtained: 

*(r. = oo) = {j /g ^ = V6x(odd,odd) (?g) 

Between the two limits, the occupation numbers Pk of 
the two identified PMWMs (r s = 10 and r s = 40) have a 
mixed character, where both the "solid" and the "liquid" 
patterns are visible when one uses a logarithmic scale. This 
is what one should expect when a "supersolid" is forming: 
the usual spreading of the "Fermi sea" being accompanied 
by the gradual emergence of the Fourier spectrum of the 
"Wigner solid" . 
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